The analysis of single crystal diffuse scattering using the Reverse Monte Carlo 

method: Advantages and problems 
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The scattering from crystals can be divided into two parts: Bragg scattering and diffuse scattering. 
The analysis of Bragg diffraction data gives only information about the average structure of the 
crystal. The interpretation of diffuse scattering is in general a more difficult task. A recent approach 
of analysing diffuse scattering is based on the Reverse Monte Carlo (RMC) technique. This method 
minimises the difference between observed and calculated diffuse scattering and leads to one real 
space structure consistent with the observed diffuse scattering. 

The first example given in this paper demonstrates the viability of the RMC methods by refining 
diffuse scattering data from simulated structures showing known occupational and displacement 
disorder. As a second example, results of RMC refinements of the diffuse neutron- and X-ray 
scattering of stabilised zirconia (CSZ) are presented. Finally a discussion of the RMC method and 
an outlook on further developments of this method is given. 
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I. INTRODUCTION 

Crystal structure analysis based on Bragg diffraction 
data reveals only information about the average crystal 
structure, such as atomic positions, thermal ellipsoids 
and site occupancies. Any departure from the strictly 
long-range-ordered average structure of the crystal gives 
rise to diffuse scattering containing information about 
static or thermal disorder within the studied material. 
Since many properties like optical properties, hardness, 
ionic conductivity are governed by structural disorder, 
the determination of the defect structure based on diffuse 
scattering data can reveal important information about 
the studied material. The development of area detec- 
tors (proportional counters, image plates, CCD) in recent 
years has enormously increased the capability for measur- 
ing generally weak diffuse neutron and X-ray scattering 
data. However, the interpretation and analysis of diffuse 
scattering remains a generally difficult task. The avail- 
ability of modern (super)computers opens a wide area 
of computer simulation techniques to aid the analysis of 
diffuse scattering. An overview of traditional approaches 
and the analysis of diffuse scattering via computer sim- 
ulations can be found in JlJ, further general information 
about disorder diffuse scattering can be found in numer- 
ous review articles |@-0|. We want to focus in this paper 
on the analysis of diffuse single crystal scattering using 
the RMC method developed by McGreevy and Pusztai 

The RMC refinement technique minimises the differ- 
ence between observed and calculated diffuse scattering 
intensities as a function of the positions and occupancies 
of the atomic sites in the model crystal. Although the 
RMC method is known for about 10 years, the applica- 
tion of RMC to diffuse single crystal scattering data was 
reported just recently in a neutron diffraction study of 
the diffuse scattering of ice Ih [p|. The advantage of the 
RMC method is the fact, that it is a model free method 



to analyse diffuse scattering, i.e. no assumptions about 
the particular disordered structure under investigation 
have to be made. In general a RMC simulation gives 
one real-space structure consistent with the experimen- 
tal data. The remaining difficulty is the interpretation of 
the resulting structure. 

The aim of this paper is to give an introduction into 
the RMC simulation technique for the analysis of diffuse 
single crystal scattering and to discuss the advantages 
and difficulties of the RMC method. 



II. THE RMC METHOD 

In general Monte Carlo methods can be described as 
statistical simulation methods involving sequences of ran- 
dom numbers to perform the simulation. In the past sev- 
eral decades this simulation technique based on the algo- 
rithm developed by Metropolis [ 10| has been used to solve 
complex problems in nuclear physics, quantum physics, 
chemistry as well as for simulations of e.g. traffic flow or 
econometrics. The name Monte Carlo was coined during 
the Manhattan Project in World War II, because of the 
similarity of statistical simulation to games of chance, 
and because the capital of Monaco was a center of gam- 
bling. In this analogy the 'game' is a physical system 
and the scientist might 'win' a solution for his partic- 
ular problem. An excellent application for this kind of 
statistical method is the study of diffuse scattering and 
subsequently the solution of the underlying defect struc- 
ture. One possible approach is the (direct) Monte Carlo 
(MC) modeling of a defect structure from a given set of 
near-neighbour interaction energies |IJ. The same basic 
algorithm is used for RMC simulations to minimise the 
difference between calculated and measured diffuse scat- 
tering as described in the following section. 
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A. How does it work ? 

As described in the introduction, the aim of the RMC 
simulation process is to minimise the difference between 
observed and calculated diffraction pattern. As a first 
step, the scattered intensity is calculated from the cho- 
sen crystal starting configuration and a goodness-of-fit 
parameter x 2 is computed. 
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The sum is over all measured data points hj, I e stands 
for the experimental and I c for the calculated intensity. 
The RMC simulation proceeds with the selection of a 
random site within the crystal. The system variables 
associated with this site, such as occupancy or displace- 
ment, are changed by a random amount, and then the 
scattered intensity and the goodness-of-fit parameter x 2 
are recalculated. The change A% 2 of the goodness-of-fit 
X 2 before and after the generated move is computed. Ev- 
ery move which improves the fit (A^ 2 < 0) is accepted. 
'Bad' moves worsening the agreement between observed 
and calculated intensity are accepted with a probability 
of P = cxp(— A% 2 /2). As the value of A% 2 is propor- 
tional to 1/cr 2 , the value of a has an influence on the 
amount of 'bad' moves which will be accepted. Obvi- 
ously there are two extremes: For very large values of 
cr, the experimental data are ignored (x 2 ~ 0) and with 
very small values of a the fit ends up in the local min- 
imum closest to the starting point, because there is a 
negligible probability for 'bad' moves. The parameter a 
acts like the temperature T in 'normal' MC simulations. 
The RMC process is repeated until x 2 converges to its 
minimum. 

The result of a successful RMC refinement is one real 
space structure which is consistent with the observed dif- 
fuse scattering data. In order to exclude chemically im- 
plausible resulting structures additional constrains, e.g. 
minimal allowed distances between atoms, may be intro- 
duced. 



B. RMC software 

The program used for the RMC refinements presented 
in this paper is DISCUS [pLlfl - For practical use it is nec- 
essary to include a scaling factor / and a background 
parameter b in the previous definition of the goodness- 
of-fit x 2 ■ A weight w(h) is included as well. DISCUS 
allows the user to choose a particular weighting scheme 
or read weights from a separate input file. The definition 
of x 2 becomes: 



W(hi)[l e (hi) 

1=1 



(f-iM + b)}' 



(2) 



o • 

\ 




• 


• 


• • 


► 


• 


o 


(b) shift atom 
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(c) swap displacements 
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FIG. 1. Illustration of the different RMC modes: (a) swap 
atoms, (b) shift atom and (c) swap displacements from aver- 
age position. 

As before I e (hi) stands for the measured intensity at 
the reciprocal point h^, and J c (hi) is the calculated in- 
tensity in that point. The summation is over all N ex- 
perimental data points. Three different ways to calculate 
the scale / and background b are implemented. First the 
user can define fixed values for both: f — fa,b — bg. Sec- 
ondly, the background can be set to a fixed value b = bo 
and the scaling factor / is computed according to: 
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Alternatively both values / and b can be refined dur- 
ing the RMC refinement. Equation (||) shows the corre- 
sponding definitions. 

The parameters / and b are computed in each RMC 
cycle and have usually large starting values as long as 
there are big differences between calculated and observed 
data. After every RMC move the resulting scattering 
intensity and the x 2 value is calculated. In order to 
save computing time only the contribution of the mod- 
ified atoms to the scattering is calculated. The differ- 
ence Ax 2 = Xoid — Xnew 1S taken to decide if the move 
will be accepted or not as described in the previous sec- 
tion. The program calculates separate scaling factors and 
background parameters for every used plane of experi- 
mental data. 

The program DISCUS is capable of modeling occupa- 
tional as well as displacement disorder and so far we have 
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called each crystal modification simply RMC move. In 
practice we use a mode ('switch-atoms') of simulation 
in which occupational disorder is modeled by swapping 
two different randomly selected atoms (Fig. Q). This 
procedure forces the relative abundances of the different 
atoms within the crystal to be constant. It should be 
noted, that vacancies are treated as an additional atom 
type within the program DISCUS. The introduction of 
displacement disorder is realised in two different ways. In 
the first method a randomly selected atom is displaced by 
a random Gaussian distributed amount ('shift', Fig. |l|b). 
Alternatively the displacement variables associated with 
two different randomly selected atoms are interchanged 
('switch-displacements', Fig. [l]c). The latter method has 
the advantage that the overall mean-square displacement 
averages for each atom site can be introduced into the 
starting model and these will remain constant through- 
out the simulation. Additional information about the 
program DISCUS can be found on the World- Wide- Web 



III. EXAMPLES 
A. RMC test refinements 

In this first example, we test the viability of the RMC 
simulation technique for systems showing occupational 
and displacement disorder in combination. First a dis- 
ordered structure with given correlation parameters and 
displacements is created and the diffuse intensity calcu- 
lated. This intensity is used as input for the RMC refine- 
ment. Finally the resulting structure is compared to the 
expected disordered structure. Details about the com- 
plete series of test simulations can be found in (l2) . 

The simulated structure used to calculate the 'experi- 
mental' data for the RMC refinements was a 2D square- 
symmetric crystal with a size of 50x50 unit cells, one 
atom (Zr) at (0,0,0) with an occupancy of 0.83 and a lat- 
tice constant of a=5A. We chose the defect 'test' struc- 
ture to consist of preferred vacancy pairs in <11> di- 
rection and a subsequent relaxation of the surrounding 
atoms towards the vacancy. The vacancy ordering can 



Input 



Run A 



Run B 



Run C 



ClO 
Cll 

dzr — zr 
^r-uac[-^] 

R 



-0.203 
0.523 
5.05(12) 
4.89(9) 



0.188 
0.290 
5.02(9) 
4.95(11) 
39.3 % 



-0.076 
0.412 
5.02(8) 
4.97(10) 

39.2 % 



-0.124 
0.465 
5.03(6) 
4.94(10) 
7.7 % 



TABLE I. Input structure and results of RMC test refine- 
ments 



be described using the correlation coefficient Cjj which is 
defined as: 
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Pij is the joint probability that both sites i and j are 
occupied by the same atom type and 9 is its overall oc- 
cupancy. Negative values of Cy correspond to situations 
where the two sites i and j tend to be occupied by differ- 
ent atom types while positive values indicate that sites 
i and j tend to be occupied by the same atom type. A 
correlation value of zero describes a random distribution. 
The maximum negative value of Cy for a given concen- 
tration 9 is —9/(1 — 9) (Pij — 0), the maximum positive 
value is +1 (Py = 9). We will refer to the correlation 
coefficient of nearest neighbours in <10> direction as cio 
and in <11> direction as c\\. The achieved vacancy con- 
centration 9 and correlation values cio, en are listed in 
Table |. The resulting structure is characterised by a pos- 
itive correlation en and a negative value for cio close to 
its maximum negative value of —0/(1 — 9) = —0.205. In 
other words, the vacancy ordering is given by preferred 
<11> vacancy pairs and avoided <10> pairs. 

The simulated structure and the corresponding diffrac- 
tion pattern are shown in Figure |a and |d, respec- 
tively The scattering pattern was computed on a grid of 
301x301 points for neutron scattering at a wavelength of 

lA. 

The RMC simulations started from a structure show- 
ing a random vacancy distribution and displacements 
corresponding to an isotropic temperature factor of 
B=0.5 A 2 . Simulation A was carried out by alternately 
executing one cycle in 'swap atoms' mode followed by one 
cycle in 'swap displacements' mode. A total of 15 cycles 
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FIG. 2. Input structure (a) for RMC test refinements showing preferred <11> vacancy pairs and corresponding diffraction 
pattern (d). Resulting defect structure run A (b) and run B (c) and corresponding diffraction patterns (e) and (f). 



of each RMC mode was computed. In this paper a cycle 
will be the number of RMC moves necessary to visit each 
crystal site once on average. The resulting correlations, 
displacements and R- values are listed in Table |. The re- 
sulting final distortions are consistent with the expected 
size-effect like displacements, i.e. the atoms surrounding 
a vacancy are shifted towards the vacancy. The calcu- 
lated diffraction pattern (Fig. ||e) shows a satisfactory 
agreement with the experimental data at high Q values, 
where as diffuse features at positions (|, |), with h and 
k integer, at low Q are barely reproduced. These diffuse 
intensities are mainly caused by the vacancy ordering. 
The correlation values achieved (Tab. [jj) show positive 
values in the <10> and <11> directions rather than the 
expected negative correlation in <10> direction. Conse- 
quently the resulting structure (Fig. ||b) shows a clus- 
tering of the vacancies rather than the expected vacancy 
ordering present in the model structure. It appears that 
the dominant part of the diffuse scattering caused by the 
displacements has a too strong an influence on the cor- 
relations that are achieved. 

In order to model both parts of the defect structure 
more simultaneously, during run B the 'swap atoms' and 
'swap displacements' modes were alternated every 0.1 cy- 
cles, i.e. the number of moves necessary to visit 10% of 
all atom sites on average. Additionally, the data set used 
for the 'swap atoms' moves was limited to a range of 
Q| < 4.0 r.l.u., i.e. only the low angle part, less af- 
fected by diffuse scattering due to distortions, was used 
to refine the vacancy ordering. The final values for the 



distortions (Tab. ||) are similar to those of the previous 
run. However, visual inspection of the calculated diffrac- 
tion patterns (Fig. ^f) shows a much better agreement, 
especially at low Q values, with the input data set com- 
pared to run A, although the R-value has only slightly 
improved. The resulting structure (Fig. ||c) and corre- 
lation values achieved (Tab. reflect the significantly 
improved description of the simulated disordered struc- 
ture. 

The simulation results presented here indicate that the 
RMC simulation technique is a powerful tool to analyse 
diffuse scattering and to obtain information about the 
defect structure even of quite complex systems showing 
both occupational and displacement disorder. However, 
the resulting correlation and displacement values are still 
significantly smaller than the expected values present in 
the input structure. Current efforts are to improve the 
results of RMC refinements by using a different way to 
calculate the diffuse diffraction pattern. Run C in Table | 
shows the results achieved using this modified technique. 
More details can be found in [fl4f| and in the discussion 



of the RMC method in section [V A 



B. Stabilised zirconia 

The cubic phase of pure ZrOi is thermodynamically 
stable only at temperatures above 2370C. This phase, 
however, can be stabilised at room temperature by dop- 
ing with oxides of a variety of lower valent metals, e.g. 
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CaO, MgO, Y2O3. The average structure is of flue-rite 
type, space group Fm3m, with zirconium on (0, 0, 0) and 
oxygen on (4, 4). The dopant cation occupies a zir- 
conium site and in order to maintain charge neutral- 
ity, a corresponding number of oxygen vacancies is in- 
troduced resulting in important ceramic and ionic con- 
duction properties of these materials. The diffuse scat- 
tering of these materials has been investigated by numer- 
ous workers (see references in [[l5]jl|] . More recently, the 
authors of (l7| described the diffuse X-ray and neutron 
scattering of Ca-CSZ (Zro.85Cao.15O1.g5) by a model of 
correlated microdomains using a formula for the diffuse 
intensity given by jl8|,|l9|]. The model consists of two 
types of microdomains, one based on a single vacancy 
with relaxed neighbouring atoms, the other based on a 
pair of vacancies separated by | <111> over those oxy- 
gen cubes containing a metal atom. A quite different 
'modulation wave' approach to model the vacancy distri- 
bution by [po| provided further evidence for the existence 
of § <111> vacancy pairs. A model for the diffuse scat- 
tering of yttrium stabilised zirconia (Y-CSZ) with a com- 
position of Zr . 6iio. 390i. 805 based on Monte Carlo (MC) 
simulations was given by |l6|pl| ] Our most recent work 
on the diffuse scattering of CSZs p2l uses the RMC sim- 
ulation technique to analyse the complex defect structure 
of these materials. 

These RMC simulations of Ca-CSZ used neutron as 
well as X-ray diffuse scattering data as input. The X-ray 
data were collected on our PSD diffractometer system 
p3| . The neutron diffraction data used in this study were 
collected by |l8| . A total of three layers of X-ray diffuse 
scattering, i.e. 0.3c*, 0.5c* and 0.7c*, and two layers of 
neutron diffuse scattering, i.e. 0th and 2nd layer of [110]- 
zone, were used as input for the RMC simulations (two 
layers are shown in Fig. ^a,b). Furthermore, the cubic 
symmetry of the crystal was taken into account assum- 
ing, that all symmetrically equivalent planes contain the 
same data as the experimental planes. This resulted in 
an effective number of data points in excess of 900000. 
The first RMC simulations were carried out using the 
largest computationally feasible crystal size of 20x20x20 
unit cells containing a total of 96000 atoms (including va- 
cancies). A total of 5 RMC cycles was carried out, giving 
a R- value of 29.6% (Run A). However, the resulting cor- 
relations (Tab. ||) suggest that the good fit was mainly 
obtained by longer ranging correlations rather than the 
local disorder we are interested in. 

In order to force the RMC process to model the diffuse 
scattering using correlations on a local scale, a series of 
RMC refinements, presented here, was carried out based 
on a crystal of only 5x5x5 unit cells in size. A total of 
6 RMC refinements were computed (Run B), each run 
iterated for 25 cycles. The oxygen-oxygen vacancy or- 
dering and the Zr-Ca ordering as well as displacements 
for all atoms were modeled. The simulation mode was 
similar to the one used for the test simulations described 
in the last section. The resulting average R-value for all 
refinements is 33.1%, the average of the resulting diffuse 



Neighbour Run A Run B 

Correlations 



VAC 


- VAC : ±(100) 

- VAC : ±(110) 

- VAC : |{111)* 


-0.005 


-0.008(11) 


VAC 


-0.014 


-0.011(10) 


VAC 


0.006 


-0.009(17) 


VAC 


- VAC : |(111) 


0.008 


0.015(7) 


Ca- 


Ca : |(110) 


-0.037 


-0.052(33) 


Ca- 


Ca : (100) 


-0.009 


-0.008(17) 



Displacements [A] 



VAC-O 


: 4(100) 


-0.011 


-0.031(13) 


VAC-O 


: f (110) 


0.003 


0.009(5) 


VAC-O 


: 1(111) 
±(110) VAC 


0.005 


0.014(7) 


Zr- Zr: 


0.024 


0.064(39) 


Zr - Zr : 


1(110)° 


-0.002 


-0.012(6) 


Zr-Ca: 


±(110) VAC 
1(110)° 


0.020 


0.029(35) 


Zr-Ca: 


-0.007 


-0.019(10) 



TABLE II. Results of the RMC refinements of the diffuse 
scattering of Ca-CSZ 



neutron and X-ray patterns of all RMC refinements are 
shown in Figure |^c and d, respectively. 

Two different types of correlations are listed in Ta- 
ble |l|: The oxygen vacancy- vacancy correlations are rep- 
resented by 'VAC- VAC and given for nearest neigh- 
bours (| <100>), next-nearest neighbours (1 <110>) 
and second-nearest neighbours over those oxygen cubes 
containing no cation (| <111>*) and those cubes filled 
with a cation (| <111>). Additionally the Ca-Ca occu- 
pancy correlations for nearest neighbours (5 <110>) and 
next-nearest neighbours (<100>) arc given. Inspection 
of Table |l| shows that the average of all RMC refine- 
ments using the 5x5x5 unit cell crystal (Run B) results 
in an oxygen- vacancy ordering scheme similar other mod- 
els proposed in the literature. However, the given stan- 
dard deviations indicate there is a large variation in the 
values obtained from the different RMC refinements, so 
that the actual values are barely significantly different 
from zero. The RMC run A using the 20x20x20 unit 
cell crystal shows the same general trends in the cor- 
relation values, although here magnitudes are generally 
even lower and in one case of opposite sign (VAC- VAC 
1 <111>*). Despite the actual values being small, the 
trends are quite definite and reproducible. It should be 
noted, that the larger values are obtained using a smaller 
crystal size, which is consistent with the view that the 
smaller system is less able to use long-range correlations 
to obtain the fit. The Ca-Ca correlations show negative 
values for nearest and next-nearest neighbours for both 
refinements, although again these are very low in magni- 
tude. Such negative correlations indicate Ca-Ca nearest 
and next-nearest neighbours tend to be avoided, i.e. less 
probable than in a random cation distribution. Again the 
smaller model crystal size leads to slightly larger negative 
values. 

The displacements listed in Table || are the distances 
from the average fluorite position. The results show 
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that the nearest neighbour oxygens are shifted towards 
the vacancy along the <100>-direction whereas next- 
nearest and second-nearest neighbouring oxygen atoms 
are moved away from the vacancy. The metal-metal 
distance along | <110> is shorter than in the average 
structure if both bridging oxygen sites are occupied. If 
one or both of these bridging sites is vacant, the met- 
als are shifted further apart. The results listed in Table 
show no significant difference for the displacements of 
Zr-Zr pairs compared to Zr-Ca pairs. As the correla- 
tion values, the refinement B shows the more significant 
values compared run A using the larger model crystal 
size. However, the large standard deviation of these dis- 
placements indicates large statistical errors due to the 
small crystal size. The resulting defect structure is con- 
sistent with previously reported models of the disorder 
in CSZ materials. These calculations as well as as well as 
a study of the diffuse scattering of TlSbOGeO^ [^4| have 
shown, that the size of the model crystal is determined 
by two conflicting requirements: a large crystal size gives 
sufficiently smooth diffraction patterns and statistically 
significant correlation parameters, but fit is obtained by 
many longer-range correlations rather than by few short- 
range correlations which are the parameters we are inter- 
ested in. A small crystal size on the other hand reduces 
the number of variables and longer-range correlations but 
the calculated diffraction pattern is too noisy for a satis- 
factory fit and the resulting correlation values have large 
errors. Further discussion of this problem a nd a way to 
improve the RMC results is given in section IV A of this 
paper. 



IV. DISCUSSION 

Generally the RMC simulation technique generates a 
disordered structure consistent with the observed dif- 
fuse scattering by minimising the difference between cal- 
culated and measured diffuse scattering patterns. The 
RMC test simulations described here have demonstrated 
that the RMC technique is a viable tool to analyse dif- 
fuse scattering and to obtain information about the de- 
fect structure even for quite complex systems showing 
occupational and displacement disorder. The study of 
single crystal diffuse scattering of Ca-CSZ confirms that 
the RMC method is able to determine the characteris- 
tic features of the complex defect structure of a 'real' 
disordered system. The combination of neutron and X- 
ray diffuse scattering data allowed the analysis to include 
oxygen-oxygen vacancy ordering as well as cation order- 
ing and both oxygen and metal displacements. The re- 
sulting defect structure shows features consistent with 
results previously reported in the literature. 

However, the resulting correlation parameters and dis- 
placements for the simulations of the diffuse scattering 
of Ca-CSZ were barely significantly different from zero 
although the trend towards the particular defect struc- 



ture was reproducible. One important problem of the 
RMC simulation technique for single crystals is the size 
of the model crystal, determined by two conflicting re- 
quirements. This 'crystal size' problem will be discussed 
in more detail in the next section. 



A. The 'crystal size problem' 

RMC simulations of the diffuse scattering of CSZ pre- 
sented here as well as a study of the diffuse scattering of 
TISbOGeOi |Q have shown, that the size of the model 
crystal is determined by two conflicting requirements: a 
large crystal size gives sufficiently smooth diffraction pat- 
terns and statistically significant correlation parameters, 
but the fit is obtained by many longer-range correlations 
rather than by few short-range correlations which are the 
parameters we are interested in. A small crystal size 
on the other hand reduces the number of variables and 
longer-range correlations but the calculated diffraction 
pattern is too noisy for a satisfactory fit and the result- 
ing correlation values have large errors. 

One way to get around the described 'crystal size prob- 
lem' is the calculation of the diffuse scattering pattern as 
the average of intensities calculated from small regions 
('lots') within the model crystal chosen at random. This 
procedure results in a smooth diffraction pattern but, on 
the other hand, structural changes during the RMC re- 
finement are done on a local scale within a single lot. 
Our current efforts are to incorporate this method of 
calculating smooth diffraction patterns by averaging the 
diffuse intensity of small areas within the crystal in the 
RMC refinement process. Additional to the high qual- 
ity diffraction patterns obtained, this method results in 
a RMC refinement using only local correlations within 
one lot. First refinements of the diffuse scattering of 
the test structures |Q described in section III A| of this 
paper using 'lots' were carried out (Run C). Although 
the same disordered structure was used for these tests, 
the diffuse scattering used as 'experimental data' for the 
RMC refinement was calculated using 'lots' as well. Ob- 
viously the size of these 'lots' must be large enough to 
contain all significant interactions. The resulting struc- 
ture and diffraction pattern are shown in Figure |]a and 
4b. Inspection of Table | shows for Run C correlations 
values and displacements significantly closer to the ex- 
pected values present in the input structure. It should 
be noted, that the refinement using 'lots' requires sub- 
stantially more computer resources compared to 'normal' 
RMC refinements. These test simulations demonstrate 
that the new RMC simulation technique using 'lots' im- 
proves the resulting local disorder for the type of systems 
discussed in this paper. 

One might argue, that the desired restriction to a 'lo- 
cal' scale in fact limits the potential of the RMC method. 
However, for many systems a 'local disorder' model gives 
usually not only the simplest description of the particu- 
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FIG. 3. Diffuse scattering of Ca-CSZ: (a) neutron scattering of 0th layer of the [110]-zone and (b) X-ray scattering of layer 
0.5c*. Corresponding resulting diffraction patterns from RMC refinement: (c) neutron and (d) X-ray. 




FIG. 4. RMC refinement results from run C using 'lots': 
(a) resulting structure and (b) corresponding diffraction pat- 
tern. 



lar disorder, but the local chemistry and near-neighbour 
interactions are frequently the interesting properties of 
the found defect structure. The test examples (see sec- 
tion III A) showing mainly short range order, were sig- 



nificantly better refined using the modified RMC simula- 
tion method using 'lots' compared to the 'normal' RMC 
refinement. One should bear in mind, that the RMC 
method produces the most disordered structure consis- 
tent with the data [01. 



B. What next ? 

The RMC simulation technique has a large potential 
for further improvements and developments and, in the 
authors opinion, will continue to play its important role 
as method to analyse single crystal diffuse scattering. 

Further developments of the RMC simulation tech- 
nique to analyse diffuse scattering of single crystals are 
twofold: First constrains like minimal allowed distances 
between atom types can be used to avoid resulting struc- 
tures which are unlikely or even impossible from a chem- 
ical of physical point of view. Those constrains can be 
highly specialised for a particular problem. A second 
major improvement can be expected from RMC simu- 
lations combining single crystal diffuse scattering data 
with other experimental data, e.g. powder diffraction 
data or EXAFS data. Thus a large RMC model could 
be constrained using experimental data more dependent 
on short range order. An the other hand unwanted long 
range fluctuations of might be suppressed by including 
Bragg data in the RMC refinement process. Besides 
the options already mentioned, one significant improve- 
ments would certainly be the possibility to used fully 3- 
dimensional data sets of diffuse scattering for the RMC 
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refinement which is currently beyond available computer 
resources. With the continuing increase of available com- 
puting power, we will certainly see further developments 
in the area of computer simulations to aid the analysis 
of diffuse scattering in general and the RMC method in 
particular. 
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